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Abstract: We present a semi-analytic method for the integrand reduction of one-loop 
amplitudes, based on the systematic application of the Laurent expansions to the integrand- 
^ , decomposition. In the asymptotic limit, the coefficients of the master integrals are the 

j^ \ solutions of a diagonal system of equations, properly corrected by counterterms whose 

parametric form is known a priori. The Laurent expansion of the integrand is implemented 
through polynomial division. The extension of the integrand-reduction to the case of 
numerators with rank larger than the number of propagators is discussed as well. 
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1. Introduction 

The recent development of novel methods for computing one-loop scattering amplitudes 
has been highly stimulated by a deeper understanding of their multi-channel factorization 
properties under special kinematics enforced by on-shellness [1,2] and generalized unitar- 
ity [3,4]. Analyticity and unitarity of scattering amplitudes have then been strengthened 
by the complementary classification of the mathematical structures present in the residues 
at the singular points. They naturally arise after uncovering a relation between numerator 
and denominators of one-loop Feynman integrals, yielding the multipole decomposition 
of Feynman integrands originally proposed in a four-dimensional framework by Ossola, 
Papadopoulos and Pittau (OPP) [5,6], and later extended to dimensionally regulated am- 
plitudes by Ellis, Giele, Kunszt and Melnikov (EGKM) [7-9]. 

The use of unitarity-cuts and complex momenta for on-shell internal particles turned 
unitarity-based methods into very efficient tools for computing scattering amplitudes. 
These methods, recently reviewed in [10-12], exploit two general properties of scatter- 
ing amplitudes, such as analyticity and unitarity: the former granting that amplitudes can 
be reconstructed from the knowledge of their (generalized) singularity-structure; the latter 
granting that the residues at the singular points factorize into products of simpler ampli- 
tudes. Unitarity-based methods are founded on the underlying representation of scattering 
amplitudes as a linear combination of master integrals (Mi's) [13,14], and their principle 



is the extraction of the coefficients entering in such a hnear combination by matching the 
cuts of the amphtudes onto the cuts of each MI. 

The multi-particle pole decomposition for the integrands of arbitrary scattering ampli- 
tudes emerges from the combination of analyticity and unitarity with the idea of a reduction 
under the integral sign. 

The principle of an integrand-reduction method is the underlying multi-particle pole 
expansion for the integrand of any scattering amplitude, or, equivalently, a representation 
where the numerator of each Feynman integral is expressed as a combination of products 
of the corresponding denominators, with polynomial coefficients. Each residue is a (multi- 
variate) polynomial in the irreducible scalar products (ISP's) formed by the loop momenta 
and either external momenta or polarization vectors constructed out of them. The inte- 
grand reduction method has been recently shown to be applicable to scattering amplitudes 
beyond one-loop as well [15,16]. 

The polynomial structure of the multi-particle residues is a qualitative information that 
turns into a quantitative algorithm for decomposing arbitrary amplitudes in terms of Mi's 
at the integrand level. In the context of an integrand-reduction, any explicit integration 
procedure and/or any matching procedure between cuts of amplitudes and cuts of Mi's 
is replaced by polynomial fitting, which is a simpler operation. Within this algorithm 
the (dimensionally regulated) integrand of a given scattering amplitude is the only input 
needed for sampling the integrand on the solutions of generalized on-shell conditions. 

The original algorithm [5,7] is based on the solution of a triangular system of equations 
to be solved top-down, from the determination of the 5-point coefficients to the 1-point 
coefficients. At any step of the reduction, the Gauss-substitutions requires the subtraction 
of the of the coefficients determined in the previous steps. Accordingly, the algorithm 
proceeds by subtractions at the integrand level, requiring the knowledge of the already 
reconstructed residues. In other words the determination of the coefficients of the n-point 
residues (1 < n < 5), requires the subtraction of the residues of m-point functions with 
n < m < 5. The integrand reduction method has been implemented in the publicly 
available libraries CutTools [17] and samurai [18], as well as in several multi-purpose 
codes [19-30]. 

In this paper, we exploit the asymptotic behaviour of the integrand decomposition, and 
propose a simpler technique for the integrand-reduction of one-loop scattering amplitudes in 
dimensional regularization. The use of Laurent series within the unitarity-based methods 
has been mainly developed in the context of analytic calculations [31-35]. 

Elaborating on the the techniques proposed by Forde [33] and Badger [35] we apply sys- 
tematically the series expansion to the integrand decomposition formula of OPP/EGKM. 
We show how the advantages of the analytic techniques can be incorporated in a refined 
semi-analytic algorithm which determines the coefficients using Laurent expansions rather 
than polynomial fitting. The main features of this algorithm are the following. 

• The coefficients of 5-point functions are never needed and do not have to be computed. 

• The spurious coefficients of the 4-point functions do not enter the reduction and are 
not computed. The rational terms coming from higher-dimension 4-point functions 



can be computed analytically from quadruple-cuts only. 

• The computation of 3-, 2-, and 1-point coefficients is independent of the residues of 
the 4-point functions. In particular the 3-point coefficients are computed from triple 
cuts, without any subtraction. 

• The subtraction at the integrand level is replaced by the subtraction at the coefficient 
level. Indeed in the original reduction method the subtractions guarantee the poly- 
nomiality of each residue allowing its determination through polynomial fitting. The 
Laurent expansion makes each function entering the subtraction separately polyno- 
mial. Therefore the subtraction of higher point residues can be omitted during the 
reduction. Its effect is accounted for correcting the reconstructed coefficients. 

• The correction terms of 2-, and 1-point functions are parametrized by universal func- 
tions in terms of the higher points coefficients. 

• The application of the Laurent expansion for the determination of 3-, 2-, and 1-point 
coefficients is implemented via polynomial division. The Laurent series is obtained 
as the quotient of the division between the numerator and the product of the uncut 
denominators, neglecting the remainder. 

This algorithm has been implemented in CH — \- and in MATHEMATICA using the S@M 
package [36]. The semi-analytic implementation in CH — \- has been designed as a reduction 
library to be linked to codes like GoSam and FormCalc which provide analytic expression 
of the integrands, or to any package which can provide the tensor-structure of the integrand 
as described in [37,38] and in the more recent Open-Loop technique [39]. The version in 
MATHEMATICA has been used to obtain closed formulas for the coefficients, depending 
on the vector-basis associated to each cut and on the generic tensors appearing in the 
integrand. In this latter fashion, the reduction procedure is replaced by a simple pattern- 
matching. 

2. Integrand decomposition 

In this section we collect the relevant formulae of the integrand reduction methods intro- 
duced in [5-9] following the notation of [18]. 

The reduction method is based on the general decomposition for the integrand of a 
generic one-loop amplitude. Any one-loop n-point amplitude can be written as 



An= d'^q A{q, e) ., 

M{q,e) 



UqL>i ■ ■ ■ Un~l 
D^ = {q+p,f-m^,={q+p^f-mj-fi^ {po ^ 0) . (2.1) 

We use a bar to denote objects living in d = 4 — 2e dimensions, following the prescription 

i=i + ii , with f = q^-fj.^ . (2.2) 



The most general numerator of one-loop amplitudes M{q, e) can be thought as composed 
of three terms, 

N{q, e) = No{q, f?) + eNi{q, ^2) + e^N2{q, fi^). (2.3) 

The coefficients of this e-expansion, A'^o, A'^i and A'^2, are functions of q'^ and /x^. In the 
following discussion we will denote by A^ any element of the set {A'o, A^i, A'^2}- 

The numerator N(q,fj,'^) can be expressed in terms of denominators Di, as follows 

n— 1 n—1 n—1 n— 1 

N{q,^?)= ^ /\ijMm{q,l^^) n Dh+^/^^jM{q,^^^) W Dh+ 

i«m hyti,j,k,i,m i«£ h^i,j,k,i 

n—1 n—1 n—1 n—1 n—1 n—1 

+ Y,/\ijk{q,l^^) n Dh + Y,^^J{q,^l^) n Dh + Y,l^i{q,^l^)\{Dh ,{2A) 

i«k hy^i,j,k i<j h^i,j i h^i 

where i « m is the lexicographic ordering i<j<k<£<m. The functions A{q,fj?) 
are polynomials in the components of q and in fi^. The decomposition ( |2.4| ) expose the 
multi-pole nature of the integrand 

n—1 A f 2\ n—1 A / 2\ n—1 \ / 2\ 

^gA^^gA^ (2.5) 

For each cut {ijk • • • ), obtained setting Di = Dj = D^ = • • • = 0, we introduce a basis 
of four massless vectors 

giijk-) ^ |ep-),eP-\eP-\eP-)} , (2.6) 

such that 

;• M =0, el -63 ' = el -el =0, 

Jijk-) (ijk-) _ (ijk-) (ijk---) _ ^ Jijk-) Jijk-)_ (ijk---) (ijk---) _ . 

^2 ■ ^3 ~ ^2 ' ^4 — U , e^ ■ e2 — ^Z ' ^4 ~ ^ ■ 

The massless vectors e^'' and 63 can be written as a linear combination of the two 

external legs at the edges of the propagator carrying momentum q + pi, say Ki and K2 , 
along the lines of [18]. In the case of double-cut, Ki is the momentum flowing through 
the corresponding 2-point diagram, and K2 is an arbitrary massless vector. In the case of 
single-cut both Ki and K2 are chosen as arbitrary vectors. In the case of quadruple-cut 
(ijkl) we define 

z;(^^-'=^) = [k, . ef''^) ef''^ + \k, ■ ef '^^^j ej^'^^) . (2.7) 

The momentum K^ is the third leg of the 4-point function associated to the considered 
quadruple-cut. To simplify our notation we will omit the indices of the cut (ijk • • • ) when- 
ever possible. 



The functions /S.{q,p?') are parametrized in terms of the basis ( |2.6D and of the vec- 
tors Q: 

A.,H„^(g,^^) = cg^"^^M^ (2.8) 

A,,,,!?, ^^') = A.^.,(cz, ^^) + cjf) + cl:f V^ + 4:f V^ , (2.9) 

A,,,(g, ^2) = A«,(,, ^2) + cj^o'^) + ctft? , (2.10) 

A,,(g, ^2) = A,^(,, ^2) + cg^^) + cg^gV' , (2.11) 

A.(g,M') = c% + cf^^{{q + p^) ■ ei) + cg((9 + Pi) ' 62) 

+ c%i{q+p^)-e3)+cf\{iq+p^)-ei) . (2.12) 

For later convenience, we define the reduced polynomials A as, 

AgM(^, M^) = (4f '^ + 4f '^ ^^) (Q + p.) • -± , (2.13) 

A.^.(g, M^) = (4:f ^ + cjf V^) (^Z + P^) ■ es + (4f ) + 4y V^) (? + P.) • e4 

+ 4'f(('?+P.)-e3)' + 4'f((9 + P«)-e4)^ 

+ 4'f (('Z + P^) • es)' + 4'f ((^ + P«) • e^f > (2-14) 

A«(g,^2) ^^fe)(^ + ^^) .,2 + ,fe)((^ + p^) .,2)2 

+ 4^3^ ('Z + Pi) • 63 + cgi^ (((7 + P^) • esf 
+ c%^ (q + Pi) ■ £4 + c%\iq + p,) ■ e^f 
+ c%\iq + Pi) ■ e2){{q+p.{) ■ 63) + 4S{{q+P^) ■ e2){{q + p^) ■ 64) . (2.15) 

Neglecting terms of ©(e), the one loop amplitude can be written in terms of master 
integrals and of the coefficients of Aijkim, ^ijki, ^ijk, A^j, and Aj, 

n-l 



4 _ \^ J3M)t ,SiM)T r,,4i\ 

i<j<k<i ^ J 

n-1 . V 

_l_ V^ J i'i'JI') T _i_ (^i'i:) r [ 2i I 

i<j<k ^ J 

n-1 . ^. 

+ E |4!o^^^J + 4!lI^j[iq+P^) ■ 62] + C^^p,j[{{q+p.,) ■ e^)^] + 4^9^ ^^J [/^^'] | 



where 



n-1 

+ E4'U> (2-16) 



/,,..,Ja] ^ / A- — , /,,..,, ^ I,,..,J1]. (2.17) 



As already noted in [5,8, 12], some of the terms appearing in the integrand decomposi- 



tion (2.5) vanish upon integration. They are cahed spurious and do not contribute to the 



amphtude ( 2.16 ). Beside the scalar boxes, triangles, bubbles and tadpoles, the other master 
integrals are the linear and quadratic two-points functions [40, 41] and the integrals con- 
taining powers of fj? in the numerator. The latter can be traded with higher dimensional 
integrals [40,41] 



/„..... [(.Y/(,, .=)] = ^ n (- 3 + / <*-N-j^ 

As already noticed in [42], eq. (l2.16| ) is free of scalar pentagons. 



(2.18) 



3. Reduction algorithm 

The one loop amplitude is completely known, provided the coefficients c appearing the 
r.h.s. of eq. ( 2.16| ) are known. In this subsection we show how to get the coefficients of 
each polynomial A performing suitable series expansions. 



Quintuple Cut 

The coefficient d 



(ijklm) 
5,0 ' 



eq. (|2.8D , can be computed using a quintuple cut. However its 



actual value it is not relevant in our reduction algorithm. Therefore its computation is 
omitted. 



Quadruple cut 

The solutions of the quadruple cut D^ 



Di = 0, can be expressed as 



Q± 



(ijki) 



-pi + Xiei -I- ^262 -I- XyV ±uv±_ 



u 



'a± + -9- 



(3.1) 



where the coefficients xi,X2,Xt, and aj_ are fully determined by the cut-conditions. The 
only two coefficients that are needed are obtained from [4,35]. 



iV(g5^'''\0) 



Jijfc^) 



+ 



^(^fo-fc^)^O) 



(ijki) 



n.^,,,M D'^ (<' ^ ' 0) n.^,,,M ^^ (^-' ' ' 0) 



^{ijki) 
-4,0 ' 



N± 



Y[h^i,j,k/^h,d 



4f 'V' + oif,') . 



/x^— >oo 



Here and in the following we use the abbreviation 



(3.2) 
(3.3) 



(3.4) 



omitting whenever possible the indices of the cuts {i,j, ■ ■ ■} as well as the /i^ dependence. 



Triple cut 

The solutions of the triple-cut, Di = Dj = D^ = can be parametrized as, 

2 
(ijk) CtQ + U 

1+ = -Pi + a^iei + ^262 + tes H 64 , 

2 

(ijk) Oio + /i 

gl = -Pi + xiei + X2e2 H 63 + ^64 , 



(3.5) 
(3.6) 



where xi, X2 and ag ^-^e frozen by the triple-cut conditions. The coefficients can be obtained 
in the large-t limit according to 



N^ 



T[h^i,j,k ^h, 



'3,4 



+4TV)t+4:ft2-4:f t3+o(i). 



n 

1 
2 



hy^i,j,k 



Dh,- 



t— >oo 



t^tOO 



Sijk) Jijk) 2 



{ijk) .2 (ijk) ,3 



4,1^ + 4XV^ i + 4^2^ i^ - 4^3^ t' + o(i) , 



(3.7) 
(3.8) 



N, 



n 



/i7^ij,fc -^/i,+ 



+ 



N^ 



n 



hjti,j,k 



Dk 



t->oo 



4^f +4fV)+o(i)+l^(t). (3.9) 



The Bachmann-Landau symbol Vt{t) denotes terms which are non- negligible with respect to 
t as i — > cxD. Eq. ( |3.9D has been introduced in [33,35]. As explained in [12], the average over 
the two solutions cancel the contributions of the spurious coefficient of the boxes. Eqs. ( p.7| ) 
and ( p.^ ) determine the spurious coefficients of the triangles. They are independent of the 
spurious coefficients of the boxes which are 0(1) in the t ^ 00 expansion. 



Double cut 

The solutions of the double cut Di = Dj = are parametrized as follows 

{ij) , , / , N , . , Ai^ + /3o + Axi + /32xf 
g+ = -Pi + xiei + (ao + xiai)e2 + te^ -\ — 64 

(ij) , , r , X , ^^ + /3o + l^lXl + (^2X1 

q_ = -Pi + x\e\ -F (ao + xxa\)e2 ^ — 63 + te^ 



(3.10) 



where Oj and /3j are kinematical factors determined by the cut conditions. The coefficients 
can be extracted from the large-t expansion, 



iV4 



n-l A_R 



ri/iT^i J ^h^ 



Uhj^ij Dh,- 






A 



ijk, - 



D, 



ijk, - 



E 

k=/^i,j 



D, 



3J) I ^fo'),,2 I M) ™, I J.iJ) ™2 I 
^2,0 """ ^2,9 t^ "•" ''2,1 '''I "•" ^2,2 ■^l"r 



t— >oo 



-2,5 



t + c^/e ^ 



fo) +2 ^(jj) 



-2,8 



xit-FO 



Ju) , SiJ) ,,2 I ^fe') 



-2,0 "T "-2,9 



(«i) ^2 



+ '^2.9 ^ + '^2.1 ^1 + '^2.2 ^1+ 



t— >-oo 



Co 3 t "P Co H 6 Co n X^C ~p O* 



(3.11) 



(3.12) 



Eqs. ( ^.11 ) and ( 3.12 ) hold only if the uncut denominators are linear in t, namely 

Dh,±' ="^ 263,4 -iPh- Pi) t + 0{l) "ih^ij. (3.13) 

Therefore the momentum K2, entering in the definition of 63^4, has to be chosen so that 

(Ph - P^) ■ esA ^ (3.14) 

for all h 7^ i,j. 

The terms involving the reduced residues A^^ remove the contributions of the spuri- 
ous three-points coefficients. The treatment of the subtraction terms is thus different from 
the one proposed in [33-35], where the spurious three-points contamination is removed by 
subtracting all possible triple cuts constructed from the double cut {ij). 



We remark that in general neither 

N± 



Uh^ij Dh,± 



nor 



Dk,± 



are polynomial in t and 1/t, but only their difference so it is. Instead their Laurent 
expansion has the same polynomial structure of the r.h.s. of eqs. ( p. 11 ) and ( p. 12 ). For 
the "+" case we have 



N< 



Uh^ij Dh,+ 



(u) I (*i) 2 I (ij) , (ij) 2 , 

O2.0 ' ^2 9 /^ ' ^2 1 "^1 ' ^^2 2 "^l' 



t— ^oo 



^2,5 



Aij) +2 



{ij) 



t + a^^ii t - ^2 8 ^1* + O 



(3.15) 



ijk, 



DkA 



t— >oo 






(3.16) 



The "— " case is obtained by replacing (5,6,8) — )• (3,4,7). Therefore in our algorithm the 
coefficients a^*-'-* and b^"^^' ' can be computed separately, obtaining the coefficient c'*-'^ by 
their difference. 



-'2,m 



a. 



n-l 

(ij) _ sr^ Aij\k) 

2,m 2-^ "2,m 



(3.17) 



In other words the subtraction is implemented at the coefRcient-level rather than at the 
integrand- level. Moreover, given the known structure of A^^, the analytic expression of 
the coefficients b^'^^'^' can be computed once and for all, in terms of the 3-point spurious co- 
efficients, the corresponding basis, and the basis of the cut (ij). The actual semi-numerical 
procedure has to be applied only to the term involving the numerator, in order to determine 
the coefficients a^'^^K 



Single cut 



(3.18) 



We consider the following solution of the single cut Di = 0, 

(i) , ao + n^ 
q+ = -Pi + xxe\ H 62 

with ao fixed by the cut conditions. The coefficient c^ q is extracted from the large-xi limit, 

n-\ 



N. 



Wh+^^K 









Dj,+Dk, 



E 



D 



j,+ 



Si) 
'1,0 



O 



a;i— >oo 



+ Vt{xi) .(3.19) 



The symbol J7(xi) denotes terms which are not negligible with respect to xi as xi -^ cx3. 



Eq. ( 3.19 ) holds only if the uncut denominators are linear in xi, 

Dh^+ "^=~ 2ei • {ph - Pi) xi + 0(1) , for /i / i . (3.20) 

Therefore Ki and ii'2) entering the definition of the basis ei^2) have to be chosen accordingly. 



The contributions from the spurious two- and three-points coefficients are discarded 
subtracting the reduced residues A^ and A^^. The subtraction procedure differs from 
the one presented in [34], where the spurious 2- and 3-point contributions are removed 
subtracting the double and triple cuts constructed from the single cut (i). 



Also in this case we remark that 



A^., 



Uh^^Dh 



Dj^+Dk^^ 



and 



A^ 



Di 



are not separately polynomial in xi and 1/xi, but only their combination in Eq.( ^.19 ) so 
it is. Instead, their Laurent expansion has the same polynomial structure of the r.h.s. of 
eqs. (|3.ig| ). 



N. 



Uhyti Dh, 



ijk, 



Dj^+Dk,+ 



Xl^^OD 



Xl— >oo 



D 



'j,+ 



a% + ^ (^ ) + ^ (^i) ' 



6i:f ^ + o(j-]+n (xi) 



6iy + o(^)+f^(xi) . 



(3.21) 
(3.22) 
(3.23) 



The coefficients a^*\ JjiM^i^)^ and b^"^'^' can be computed separately, and finally the 1-point 
coefficient read, 



n— 1 n— 1 

-1,0 ~ "i,o Z^ "1,0 Z^"i,o 

j<k^i jy^i 



(3.24) 



{i) 



The semi-numerical procedure has to be used to determine the coefficient o^^ q only. Indeed, 
given the known structure of A^^, and A^, the parametric form of the coefficients b^Q 



(i\'j) 

and b^Q is universal and can be computed once and for all. 



4. Implementation 

As shown in the previous section, our method requires, from quadruple- to single-cut, 
one-dimensional asymptotic expansions. In this section we present their semi-numerical 
implementation . 

Quadruple cut 

The coefficient C4 4 is computed performing a large //^ expansion of 

(4.1) 



\.ihj^i,j,k,i^h, - 

cfr. eq. (|3.3| ). Both N and D^ are polynomial in u 



N+ = Y,hu' = N {uv^,vW) + OK-i), Dh, + = dh,o + dh,iu . (4.2 

Jijki) 

.VVV,i V^WUliJ-UiiJ-g,, UJ.1V, V,U<..-iJ.iV,iV.-lJ.U C4 4 

to fr, 

(ijkt) Jr Ji 



By power counting, the coefficient c^ 4 is non- vanishing only if r = n and it is proportional 



KY m^lj,k,i dH,i {vlY Wl~4,,kpPh ■ ^^) 



(4.3) 



The coefficient /^ can be obtained from the analytic expression of N(uv±,v'^u'^). For 
instance this procedure can be easily implemented in FormCalc and GoSam which use 
the symbolic manipulations programs mathematica and/or FORM [43]. 

Triple, double, and single cuts 

Along the reduction procedure, the integrand of the n-ple cut is a multivariate function of 
(5 — n) variables, corresponding to the parameters of the loop momentum not fixed by the 
on-shell conditions. Each expansion is performed with respect to one variable only, say r. 
The solution of the triple, of the double and of the single cut reads as follows 

'Zc^ut = %':i+< + ^<, (4.4) 

in terms of the (cut-dependent) momenta r/_i, rfQ, and r/i. For each cut the generic term 
to be expanded is a ratio of the type, 

,^'^' . (4.5) 

nt;ok(r) 

where D^ is an uncut propagators, 

2 

Dh{T) = '^^, P,(r) = V4,r\ (4.6) 

r ^ — ' 

i=0 

The function F can either be the original numerator N or any of the reduced residues A . 



10 



If F is a reduced residue A the large r expansion of eq. ( |4.5| ) is universal and can be 
performed analytically, cfr. section |3[ 
If F is the numerator A^ we have 

nr) = ^, Hr) = r'-N (g^ut, M^) = f; /^ r\ (4.7) 

^ i=0 

where r is the rank of the numerator. The coefficients fi depend implicitly on the cut 
through the momenta ??_i, ryoi and ??i, 

fi = fi{v~-i,Vo,rii)- (4.8) 

Their parametric expression can be obtained from either the tensor structure of the inte- 
grand or the analytic form of the numerators, as provided by codes like FormCalc and 
GoSam. In the F = N case the large r expansion is numerically implemented through 
polynomial division, as described below. In the following we assume r > k, otherwise the 
ratio ( [4.5| ) vanishes in the r — t- oo limit. 

Step 1. We start by dividing F by Vq obtaining 

The quotient Qo is a polynomial of degree 2r — 2, while the remainder TZq is a polynomial 
of degree one. In the large-r limit, the contribution of the latter can be neglected, since 

J^f^^ ^^~ O (^) , (4.10) 

therefore 

F{t) ,^oo Qo(r) 



rfc-l -n / N -rrfe-l . 



Step 2. We perform the division by the successive denominator 2?i 



(4.11) 



1^^ =.(.,. 5iM. ,.., 

The quotient Qi is a polynomial of degree 2r — 4 while the remainder TZi has degree one. 
In the large-r limit, the contribution of TZi drops out as well, hence 



Fir) r^oo Qi(r) 



rfc-l -T-, r \ -Tfk-l ■ 



(4.13) 



Last step. After reiterating this procedure over the remaining denominators, 1)2, ■ ■ ■ , ^fc-i, 
we get 

2(r-fe) 



F{t) 



Utl-^kir) 



Qk-i{T)= Y, ^i^'' (4-14) 
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Finally, the Laurent expansion of eq. (4^) is given by 

Y^Si+r^kr' + o(^) . (4.15) 






m-=oDH{r) r-'^nto^/^W 



j=0 



It is worth to notice that the large r expansion can be achieved by using a smaller poly- 
nomial 

(4.16) 



i=r+k 




instead of the polynomial F defined in eq. (|4.7). Indeed 




F{t) .^oo T^{r) r^ooy^ 


-kT' + Oi- 


Y\l^,D,{r) rr~^Y[t},Vn{r) U^'^'- 



(4.17) 

We have implemented this algorithm in C-|--|- and in mathematica and verified 
its correctness reconstructing the integrands of up to sixth rank 6-point functions. Two 
numerical examples are described in the Appendix ^. A complete implementation in 
GoSam and FormCalc is planned. 

5. Example: reducing a second rank 3-point integrand 

In this section we apply the reduction procedure described above considering a rank-two 
three-point integrand of the type 

N{q) _A{q-v){q-w) ^^ ^^ 



D0D1D2 D0D1D2 

where 

Do = q^-m\ Di = {q-kif-m\ D2 = {q + k2f - rr? . (5.2) 

the "external" momenta ki and k2 taken as massless. For simplicity we consider only the 
four-dimensional part of the reduction. The extension to d-dimensions is straightforward. 
For illustration purposes, we use polynomial division . 

The cut (012) 

In order to deal with relatively compact expressions we use the basis {ei, 62, 63, 64}, where 
A«_^/^ M_^M /. _ (Ii7l2] ^_ (2|7^'|1] ._. 

The basis does not fulfill the normalization conditions ei • 62 = —63 • 64 = 1. Therefore the 



formulae (3.7-^^) have to be modified performing the substitutions 



C3,i -^ (ei • e2)c3,i, ifi = 1,4,8,9; 

C3,i -^ (ei • e2)^C3,j, ifi = 2,5; 

C3,i -> (ei • 62)^03,,, if i = 3, 6 . 
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The solutions of the triple cut are 



(012) _ rn? (012) . rr? 



^+ ='''- 2{k,.k,)t ''^ ^- ='''- 2{k,-k,)t ''- ('•') 



The functions appearing the l.h.s. of eqs. (p.7l-3.9) are 



^+ = - no"" , A iM'A (2|«^11] + (2|w|l] {l\w\2]) + (llt;l2] {l\w\2] i^ + O f i ) , (5.5) 



^- = -^(^^ (^1|H2] (2|^^|1] + (2|^|1] (l|u;|2]) + (2lt;|l] (2l«;Il] t" + O (^) , (5.6) 

since all the propagators are cut. No polynomial division is needed and the coefficients can 
be immediately computed. The non vanishing ones are 

(012) _ rn^ 



2{ki ■ k2) 
(012) _ (2|t;|l](2|^|l] 
'^'^ " {k,.k,Y ' 
(012) ^ (l|t;|2](l[w[2] 

''^^'^ {ki-k^Y ■ 

The reduced residue reads as follows 



(l|w|2](2|u;|l] + (2|i;|l](l|'u;|2]), (5.7) 

(5.8) 



(5.9) 



__(012)/(2M1]\' J012)/(1M2]'' 



^U<1) = cfr ( ^ ) + 4T' ( ^^ ) . (5.10) 



The cut (21) 

The basis used for this cut is obtained from the momenta 



AT = .f + .J, K^^J^^mR, (5.n) 

and its elements are 

e^, = \{K^-K^), ei^ = l{K^ + K^), e^ = 1^1^, e^ = i^l^. (5.12) 

The element of the basis are not canonically normalized thus eqs. ( |3.11| ) and ( ^.12]) have 
to be modified performing the substitutions 

C2,j -> (ei • e2)c2,i, ifi = 1,3,5; 

C2,j -> (ei • e2)^C2,i, if i = 2, 4, 6, 7, 8 . 

The solutions of the double cut (12) are 

(21) / ^ 2\ 1 

q+ = -k.2 + xi ei + (1 - xi)e2 + 1 63 + ( -— — —r + xi - x^ I -64 , (5.13) 

q^^^> = -k2 + xi ei + (1 - xi)e2 + 1 64 + f - + xi - x^ j -63 . (5.14) 

The coefficients Cg ; C21 and Cg 5 are obtained from eq. (|3.11| ). The large t expansion is 
obtained performing the polynomial division with respect to t, along the lines of section |^. 
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Contribution of the reduced residue - The only reduced residue entering the subtrac- 
tions is 



^012, + _ , (21|0) , , (21|0) 



(2110), 



D, 



0, 



"2,0 ' "2,1 -^1 "2,5 '' 



(5.15) 



The coefficients b^"^^'^' are universal functions of the spurious coefficients and of kine- 
matic invariants. In the case of a rank-2 three-point integrand with po = they read 
as follows 



(,(21|0) 
^2,0 



-1 



J21) 



4(er-P2)2 



(Pi + 



rrio — m. 



(21) 



i) -2ao(4 -P^ 



+ 2 



.(012) 



^21) _ {012)\ (012) 



^3,1 +200(^63 -63 j 



-3,2 



(21) 
'3 ■P2 



(21) ^(012)\ (012) 
^3 ' ^3 J ^3,2 

^(21) ^(012)\ 
'3 ' ^^3 I 



+ 



(012) (012) 

C3,l -^ C3,4; C3,2 ^ C3,5; 6^ ' ^> 6^ 



3(2|7;|l](2!i(;|l] - (l|t>|2](lji(;|2] 



4(fci • k2) 



J21|0) 
"2,1 



2(ef)-:P2)2 



(21) , (21) N f 

A ' ■p2 + aie\ ' ■ p2J [e 



(21) ^(012) \ 
■ ^3 ) 



2 e 



,(21) ^(012) 
'1 ■ ^3 



+ Qi e 



(21) ^(012)\ / (21) 



I I eg ' • P2 



J21) _ (012)\ (012) 



; 



-3,2 



+ 



(012) 
C3,2 -^ C3,5; e^ 



,(012) 



(l|t'|2](l|u;|2] - (2|t>|l](2|i(;|l] 



ki ■ k2 



(012) / (21) (012)\2 (012) / (21) (012)\2 

,(2i|o)_ 3,2 \^3 -^3 ) +C35 ye^ -64 j 



2,5 



.(21 



2(er-P2) 

(l|-t;|2](l|w|2] + (2|^;|l](2|u7|l] 
4(A;i • k2) 



(5.16) 



• Contribution of the numerator - We perform the polynomial division to compute the 
large t expansion of the ratio 



N, 



Ht) 



i— >-oo 



F^{t) 



Dn 



t Voit) t Vo{t) 



(5.17) 



The function T defined in eq. ( |4.7| ) is given by J-" = t'^N^. The reduced polynomial 
J- , eq. ( [4.16 ), is obtained from T neglecting terms of 0{t'^), 



jrR(t) = f^t'' + f^t\ 



(5.18) 
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where 



h = \({Mv\l] - (1|77|2] + (2|7;|1] - (2|t;|2] j ({l\w\l] - {l\w\2] + {2\w\l] - {2\w\2] 

h = ^ ((l|t;|l] - (l[t;|2] - (2|t;|2]) (^{l\w\l] - {l\w\2] - {2\w\2]) - (2|t>|l](2|t/;|l] 

(l|t;|l] - (l|i;|2] + (2|7;|1] - (2|i;|2]) ((l|t(;|2] + (2|i(;[l]) + (v ^ w) 

The polynomial Pq in the denominator is related to the propagator Dq, +) 

Vo{t) = t Do ,+ = do,2t^ + do,it + do,o , (5.19) 

with 



1 

+ 2 



Xl. 



do,2 = do,i = -{ki ■ k2) , dofi = (fci • ^2) 



2A;i • h 



1 ■ K2 



Xl + Xl 



(5.20) 



In the notation of eq. ( [4.15 ), the result of the polynomial division reads, 






iV_ 



+ t— >CXD 



I)n 



tVo{t) 



t—¥00 



S2t + si ,(5.21) 



where 



S2 



A 
do- 



J f-ido,2 - f4do,i 
and si = -"2 = si,o + xisi^i 

do,2 



By comparing eqs. ( |3.16| ) and ( 5.21 ) we get 



(21) 
«2,0 - ^1,0 



(21) _ ^1,1 



^2,1 



{ki ■ k2 



J21) 
^2,5 



-S2 



{ki ■ /C2 



(5.22) 



(5.23) 



while Og 2 , ^2 6 ^'^'^ ^2 8 vanish. 

Computation of the coefficients - The coefficients c^^^-' are obtained subtracting the 
coefficients b^"^^'^' to the coefficients a'^^\ according to eq. (p. 17). The non-vanishing 



ones are: 



1 



(l|t;|2] {l\w\l] + {l\w\l] (2|t;|1] - (l|u;|2] (2|7;|1] 



(21) ^ 

^''° 4(^1 . k2) y 

+ (l|u;|l] (2|7;|2] - (lju'j2] (2|t;|2] - (2|t>|2] (2|u;|l] 



(2|7;|2] (2|t(;|2] (l|t;|l] (l|t(;|l] 



+ \V -(r^ W 



.(21) 



1 



^''^ 2{ki ■ k^Y 

+ (2|7;|2](2|w|l]) +[v^w 



l\w\2] (2|w|2] - (l|u|2] (l|w|l] - (l|w|l](2|t;|l] 



.(21) 

-2,5 



-1 



4(fci • k2Y 



(l|u|2] {l\w\l\ - {l\w\l] (2|t;[1] + (l|w|2] (2|7;|1] + (ll-u^jl] (2|t;|2] 
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(l|u;|2] (2|z;|2] + (2|t;|2] {2\w\l] - ^ ' ' ^'^ ' ' ^ - ^ ' ' ^,^ ' ' ^ 



+ [v -^ w) . 

.(21) 



The remaining non- vanishing coefficient, Cj 3 , is obtained in a similar way, using eq. ( 3.12| ). 
The outcome is 

cg^ = - i(^7^ ( (1|^|1] (2|^;|1] + {l\w\2] (2|t;|l] + (l|u;|l] (2|7;|2] 
+ {l\w\2] (2|t;|2] - (2lt;|2] (2|u'|l] - (l|7;|2]|(l|u;|l] 

(2|.|2|(2|..|21 _ (l|Hl]l(lk-|l]\ ^ / , „ ^N , p,24) 



The reduced polynomial A^l is given by 

(21) (21) (21) (21) _ (21) _ (21) 

^21 = — ^ —{Q + ^2) • A:i + — ^ '—{q ■ k2) 

(21) (21) (21) , , , , (21) (21) (21) , , , , 

, 4,5 -4,3 -4,1 (i|g|2] , 4,3 -4,5 -4,1 (2|(7|i] 

+ 2 2- + 2 2- (^-'^^ 

The cut (02) 

The computation of the coefficients of the cut (02) proceeds along the lines described above. 
We define the basis using 

AT = .J, AT = .f+<l^+(^, (5.26) 



and 



e^ = i^f, eii = {K^ + K^), e^ = M^, ^^ = (fl^. (5.27) 



The solutions of the cut are: 

(02) in (02) m 

q^ =xiei+te3-^^^^^^e4, g_ = ^1 ^1 - 2(fe, . fc,) t ^3 + ^ 64 ■ (5.28) 

The non-vanishing coefficients are 

(02) _ (2|i;|2] (2|«;|2] 



2(fci . k2y 

(02) _ (2|7;|2] (2|-w|l] + (2|7;|1] (2|^;|2] + (2[t;|2] (2[u'[2] 

"'■' " 2{h-k2y 

(02) _ (1[^;[2] (2|t;[2] + (l|z;|2] {2\w\2] + (2|t;|2] (2|w|2] 

2ih-k2r 
while the reduced residue reads as follows 



(5.29) 
(5.30) 
(5.31) 



aR _ „(02) . , , / (02) (02) (02) \ _ , 

^02 — ^2 1 Q ■ Ki -\- \ C21 + C2 3 -t- C2 5 I g • K: 



^(02) ^(02) ^(02) ^(02) 

+ ''' . ''' (l|g|2] + '•' , '•' (2|g|l] (5.32) 
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The cut (2) 

We parametrize the single cut solution ( p.l8| ) in the basis (5.12), 



(2) 
1+ 



-/c2 + xiei + 



m 



1 



{2ki ■ k2) xi 



62- 



(5.33) 



The large xi expansion in eq. (|3.19|) is obtained from the large xi expansion of the sub- 
traction coefficients 



^012, 



^0, +Di, 

^21, 



aR 
^02, 



xi— ^oo 



Xl— >oo 



D: 



0, 



and from the polynomial division of the ratio 



6sr+o(^)+^(xi), 
fcS'^+o(^)+^(xi), 



6S'r+o(^)+^(xi), 



(5.34) 



Do,+Di,+ I'O^ \xi 



(5.35) 



The tadpole coefficients is 

(2) _ (2) _ , (2|01) _ , (2|1) _ , (2|0) 
•^1,0 ~ "l,0 "1,0 "1,0 "1,0 



1 



(l|w|2] (2|t>|l] - (l|u'|l](2|t;|l] 



8(A:i • A;2)2 ^ 

+ (l|ii;|l] (2|t;|2] + 3(l|iy|2] (2|i;|2] + 3(2|i;|2] (2|w|l 

,,,,,, , , (2lt;l2l (2|u;l2l (l|v|l](l|u;|l 
- (l|t;|2](l|w|l] + ^ ' ' ^^ ' ' ^ + ^ ' ' ^: ' ' ^ I +(v^w 



(5.36) 



6. Extended decomposition 

In the previous sections we assumed to deal with a renormalizable theory, where the rank 
r of the numerator can not be greater than the number of external legs n. In this section 
we extend the integrand decomposition to the case where the rank r can become larger 
than n. 

In general the residue of an m-point function is a a multivariate polynomial in /i^ 
and the ISP's characterizing the residue. Each monomial has to be irreducible and its 
maximum rank has to be at most (m -\- r — n). In the following we list the irreducible 
monomial entering each cut. For later convenience we give the decomposition of g^'^ in 
terms of the basis (12. 6D and of the vectors (|2^ 



9' 



p,v 



6162 + e^e^') - (e^e^ + e^e^) , 



(65-6^ + 6^6^) + 



+ 



V, V , 



(6.1) 
(6.2) 
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Quintuple cut, {ijkim) - The only irreducible monomial is fj,^. Indeed the residue 
of the quintuple cut does not have ISP's, thus the allowed monomials are (/^^)"- 
Moreover from eq. ( |6.lD 



= [A + m^ - p2 _ 2 (g . Pi) -2{q- ei){q • 62) + 2 (g • e^){q ■ 64)]" 

= constant terms + RSP's , (6-3) 

where the abbreviation "RSP's" means "reducible scalar products". This relation 
allows to express all the powers of ^^ in terms of a particular one, (;U^)"°. As in the 
renormalizable case we choose ao = 1 in order to decouple the contribution of the 
pentagons from the computation of the coefficients of the boxes. 

• Quadruple cut, (ijki) - The irreducible monomials are 

(a*^)" ((9 + Pi) ■ v±f with /3 = 0, 1 and a = 0, 1, 2, . . . . (6.4) 

Eq. (]^^ implies 

i{q + P^)■v±? =vl (g'-2((g+p,)-ei)((9 + K)-e2)-^^^^tg^^) 
= constant terms + terms in /x + RSP's , 
therefore the terms with /3 > 2 are reducible. 

• Triple cut, {ijk) - In this case the irreducible monomials are 

{^?r Hq + Pi) ■ 63,4)^ with a, /? = 0, 1, 2, . . . . (6.5) 

The monomials containing both 63 and 64 are reducible. Indeed from eq. ( |6.1| ) 
{{q + Pi) ■ 63) {{q + Pi) ■ 64) = constant terms + terms in ij? + RSP's . 

• Double cut, (ij) - The irreducible monomials are of the type 

(/i^)" {{q + Pi) . 63,4)'' {{q + Pi) ■ e2V with a, /?, 7 = 0, 1, 2, . . . . (6.6) 

As in the previous case, the monomials depending on both 63 and 64 are reducible. 

• Single cut, {i) - The irreducible monomials read as follows 

(/z2)° {{q + Pi) . ei,2)^ (((? + p^) ■ es)"' {{q + pi) ■ e^f with a, /3, 7, 5 = 0, 1, . . . (6.7) 

Eq. (^]^) allows to write 

{{q + Pi) ■ ei) {{q + pt) ■ 62) = ((g + Pi) ■ 63) {{q + pi) ■ 64) 

+ constant terms + terms in // + RSP's . 

Therefore the terms containing both ei and 62 do not enter the parametrization of 
the residue. 
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The residues A presented in section |2| are the most general polynomials with r < n 
satisfying these requirements. Here we show, as an example, their extension to the case 
r < n + 1. In this case, the decomposition of the numerator has to be extended as follows: 

n—l n—l n—l n—1 

i«m hj^i,j,k,£,m i«£ h^i,j,k,l 

n—l n—l n—l n—l 

i«k hj^i,j,k i<j h^i,j 

n—l n—l n—l 

+ Y. Ai{q, ^?) n Dh + A(g, ^2) J] Dh , (6.8) 

i h^i h=Q 

where the polynomials A are defined as, 

AijkimiQ, fJ- ) = '^ijkemiQ-, fJ' ) ■, 

Aijkiiq, fJ-"^) = ^ijkiiq, /w^) + 4*5 Ai^ (? + Pi) ■ v± , 

Aijkiq, fJ-"^) = \jk{q, M^) + Cg/u / + 4*10 ^^ (il + Pi) ■ es)^ 
+ 4'i? m' ((9 + P^) ■e,)' + 45? ((^ + P^) ■ es)' 
+ 4:S(.{q + Pr)-e,)\ 

Aijiq, /i^) = \jiq, ^^) + M^ (^4*10 (^ + Pi) • ^2 + c^^i (g + Pi) ■ 63 

+ 4'{k9 + Pi) • 64) + 4^3 (('^ + P*) • 62)' + 4!m(('? + Pi) ■ es)' 

+ 4!S((9 + P*) • 64)' + c^2!ll(.iq+P^) ■ e2)\{q+Pi) ■ es) 
+ 4!M{q + P^) ■ e2)\{q + p^) ■ Ci) 

+ C2i\iiq + P^) ■ e2)iiq + Pi) ■ es)^ 

+ 4*,i9((9 + Pi) • e2){{q+Pi) ■ e^f , 

+ 4S(('? +P^) • 63)' + 4'8((9 +Pi) • 64)' 

+ 4io((9 + ?'i) •ei)((g+Pi) •e3) + 4*n((9 + Pi) •ei)((g+Pi) -64) 

+ 4i2((^ + Pi) •e2)((5+Pi) •e3) + 4i3((9+K) •e2)((g+Pi) -64) 

+ 4 14 M^ + 4 15 ((9 + Pi) ■ e3)((g + Pi) ■ 64) , 
A(g,//2) = co. (6.9) 



The functions A, appearing already in the case r < n, were given in Eqs. (|2.8| )" (|2.12 ). We 



observe that the polynomial residues of 4-, 3-, 2-, and 1-point function acquire a richer 
structure, and a 0-point coefficient cq does appear. The latter coefficient is needed for 
the complete reconstruction of the integrand, but it is spurious. Indeed it multiplies a 



19 



scaleless integral, which vanish in dimensional regularization. The coefficient cq is not cut- 



constructible but it can be computed by inverting eq. (6^) in correspondence to any value 
of (g, iJ?) not annihilating any propagator. It can be shown that the 0-point coefficient 
is present in the decomposition of rank-two 1-point integrals only. In the case of higher 
point integrals this term is absent, owing to mutual cancellation between different 1-point 
polynomials Aj parametrised in terms of a common single-cut vector basis {61,62,63,64}. 

We remark that according to our new algorithm, also in the case r < n + 1 the residues 
of the 5-point functions and the spurious 4-point coefficients are not needed. Moreover, 
since the 0-point coefficient is spurious, the spurious coefficients of the 1-point functions are 
not needed as well. The other coefficients can be computed performing quadruple, triple, 
double and single cuts, using the series expansions described in section y and selecting 
the appropriate terms of the series. The new coefficients are obtained including higher 
order contributions in the expansions. The nice features of the method hereby discussed 
are not spoiled by the presence of higher rank numerators, and the coefficients of the 4- 
point functions do not affect the determination of the lower-point coefficients. We checked 
the validity of our procedure reconstructing the integrands of up to seventh rank 6-point 
functions. 

The one loop n-point amplitude is obtained upon integration of Eq. (|6.8| ). The outcome 
reads, 

An + SAn, (6.10) 



where An is given in eq. ( p. 161) and the new contribution 5 An is 

n-l 

0-A.n ^ / ^ ''3,14 ■'-ijkil^ \ 
i<j<k 

+ E USs I^Ai{<l+P^) ■ e2f] + cgi /., [/.'((^ + P,,) • 62)]} 
i<j ^ J 

+ E{^S4^«[^'] + cS5^a(('Z+Pi)-e3)((g+P.)-64)]} . (6.11) 

i 

The integral Iij[{{q + Pi) • 62)^] can be obtained from the analytic expression of the rank-3 
bubble given in section 4 of [44]. The integrals Ii[fJ.'^], Ii[{{Q + Pi) • 63) ((g + Pi) ■ 64)], 
lijil^'^iil + P'i) ■ 62)], and lijkil^^] are computed in Appendix 0. 

7. Conclusions 

In this paper we presented a procedure for the semi-analytic reduction of one-loop scattering 
amplitudes in dimensional regularization, which exploits the asymptotic behavior of the 
integrand-decomposition. 

The algorithm is based on a partial reconstruction of the numerator, where the co- 
efficients of the master integrals are determined through a simplified integrand-reduction. 
Whenever necessary, the complete integrand reconstruction can be achieved as well. 
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The analytic informations allow to avoid the computation of 5-point coefficients and 
of the spurious 4-point ones. Moreover the 4-point non-spurious coefficients do not enter 
the determination of the lower-point ones. The integrand reduction algorithm is indeed 
required only for the coefficients of 3-, 2-, and 1-point functions. 

The asymptotic expansion makes both the numerator and the subtraction terms sepa- 
rately polynomial. Therefore the subtraction of higher-point residues can be omitted during 
the reduction and replaced by coefficient-level corrective terms. The latter can be deter- 
mined a priori from the Laurent expansion of the expression of the integrand-subtraction 
terms. Therefore the actual reduction algorithm applies only to the terms involving the 
numerator, whose reconstruction is achieved by polynomial division. The coefficients of the 
3-, 2-, and 1-point functions are finally obtained as trivial combinations of the coefficients 
coming out of the polynomial division and the corrective coefficients. 

This method exploits as much as possible the known analytic structure of the integrand, 
hence it relies on the analytic structure of the numerator and its asymptotic expansions, 
used as input. It has been implemented and tested in mathematica and in C++ , using 
the polynomial division. The semi-analytic implementation in C++ has been designed 
as a reduction library to be linked to codes like GoSam and FormCalc which generate 
analytic expressions for the integrands, as well as to any package providing the tensor- 
structure of the integrand. 

We also discussed the extension of integrand-reduction methods to theories allowing 
for integrands with powers of the loop momentum larger then the number of denominators. 
We explicitly presented the extended polynomials in the case of powers larger by one unit 
than the number of denominators. The advantages of the method hereby discussed are not 
spoiled by the presence of higher rank numerators. 

We are confident that the investigation of the asymptotic regimes can ameliorate the 
integrand decomposition of scattering amplitudes beyond one-loop as well. 
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Result 

contribution a_2 
SAMURAI 
new algorithm 



-185.051790779988 +i 2.1 x lO'^^ 
-185.051790779978 - i9.0 x lO'^^ 



contribution a_i 
SAMURAI 
new algorithm 



749.007288547 
749.007288566 



i 580.971272485 
i 580.971272508 



contribution oq 
SAMURAI 
new algorithm 



-724.020439522 + i 2350.630383583 

-724.020439861 + i 2350.630383784 



coefficients evaluated 

SAMURAI 

new algorithm 



461 
386 



Table 1: Integrand reduction of Ig.e- The comparison between samurai and the new algorithm 
is shown. 



A. Numerical examples 

In this appendix we present two modest numerical applications of the hereby discussed 
algorithm. We compare the results with the ones obtained using the standard d-dimensional 
integrand reduction implemented in samurai . A comprehensive comparison between the 
two algorithms is beyond the scope of this paper. We only intend to show the potential 
benefits arising from a lighter algorithm, which requires less coefficients and which uses 
subtraction at the coefficient-level rather than at the integrand-level. 

Rank six 6-point integrand 

We consider a six-point integrand of rank-6: 



'-es 



ntoA 



A/'(9,/^')=n(9 



(A.l) 



j=i 



The momenta appearing the denominator Di are 



P'l 

M 
P2 

l-l' 

P3 
Pt 
P', 



(0,0,0,0) 

(-56.6251094805816, 0, 0, -56.6251094805816) 

(-113.250218961163,0,0,0) 

(-68.5281885958052, 33.5, 15.9, 25) 

(-48.7688869887140, 21, 31.2, 25.3) 

(-27.9148705889889, 11, 13.2, 22) 



(A.2) 
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while all the masses are assumed to be vanishing. The momenta ri are given by 

< = (1.30,5.10,0.50,0.40) 
r^ = (0.80,1.00,2.30,2.50) 
r^ = (1.90,3.20,1.77,2.11) 
r1 = (3.03,1.05,2.33,1.77) 
r^ = (3.56,5.30,3.09,2.34) 

< = (7.08, 1.98, 5.30, 4.55). (A.3) 

The integrated result is given as a series in e = (4 — (i)/2, 

^^^'"T^' I d--, Xa,a = ^ + ^ + ao + 0(e) . (A.4) 

In Table H we show the numerical values of the coefficients a^ computed with the algorithm im- 
plemented in SAMURAI and the new one. The two algorithms are in good agreement, but the new 
algorithm requires the determination of 386 coefficients, instead of the 461 required by the standard 
reduction. 

We estimate the quality of the reconstruction of the numerator using the global {N = iV)— test 
described in Section 3.4.1 of [18]. For this integrand, the reconstruction of the new algorithm is two 
digits more accurate than the one performed by samurai. 

Fermion-mass dependence of the 4-photons amplitude 



Because of the presence of products of denominators in Eq.(2.4), the numerical integrand reconstruc- 
tion may become unstable if the internal masses are larger then the kinematic invariants [18,27]. 
A simple example where such a situation may occur is the four-photon scattering in QED, 

-^ 7(fci, '11)7(^2,^2) 7(^3, /13) 7(^4, /14), (A.5) 

where ki (hi) denotes the momentum (helicity) of the corresponding particle. We define s = 
(fci + fc2)^- The leading-order process proceeds via fermionic-loop, and in the case of a single 
fermion of mass m, there are three independent helicity amplitudes 

^(+, +,+,+), Ai+, +,-,+), ^(+, +,-,-), (A.6) 

which are known analytically [45]. Mutual cancellations among contributions from different Mi's 
render the numerical evaluation of these amplitudes unstable, in particular when the ratio w? j s 
becomes large. In this Appendix we explore the behaviour of the new algorithm in this kinematic 
regime, by comparing it to the one implemented in samurai . We consider the phase space point 

fcf = (-7.0,0.0,0.0,-7.0) 

A:^ = (-7.0,0.0,0.0,7.0) 

A:^ = (7.0, 6.1126608202785198, -0.8284979592001092, 3.3089226083172685) 

fc^ = (7, -6.1126608202785278, 0.8284979592001093, -3.3089226083172703) , (A.7) 

and we vary the numerical value of m. In Figure yj we plot the relative difference (5, defined as 

A -A 



8 = 



A 



(A.i 



as a function of m} j s. .4ana (Aum) is the analytical (numerical) value of the amplitude A. For 
each helicity configuration the new algorithm seems to be less affected by this kind of inaccuracy 
than the one currently implemented in samurai . 
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Figure 1: Leading order contribution to the amplitude A. For each hehcity configuration we plot 
the quantity 8, eq. ( |A.8| ), as a function of rn^ / s. The numerical evaluation of A has been performed 
using either samurai or the new algorithm. 



B. Higher-rank integrals 



In this appendix we compute the higher-rank master integrals appearing eq. ( |6.11 ) , 



h[^l% h[{{q + P^) ■ ez){{q + p^) ■ ei% h,[^l^ {{q + p,) ■ 62% h,k[A- (B-l) 



The strategy is outlined in the Appendix I of [46]. The integrals (B.l) are obtained by projections, 
namely contracting appropriate tensors with the covariant decomposition of suitably chosen tensor 
integrals. 

For later convenience we split the d-dimensional metric tensor g^'' into the 4-dimensional part. 
g^^" , an the (— 2e)-dimensional part, g'^'^. We have 



9 — .9 +5 1 9 9ti.y — ^Cj 9ti.y1 1 ~ M i 9tJ.v'^ — ^i 



(B.2) 



for each 4-dimensional vector A:^. 
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Integral Ii[fi'^] 

The covariant decomposition of a rank-2 tadpole is 

ir ^ f d'^q^l±^^l^±^ = 5-Aoo (B.3) 

After the contraction with g^'^ we have 

/.[m^] -Jd'q^: = 2eAoo = ^ + ^(e). (B.4) 

The analytic expression of Aqo and of the other tensor coefficients appearing in this section can be 
found in [44]. 

Integral Ii[{{q + Pi) ■ es){{q + pi) ■ 64)] 



The expression of this integral can be obtained by contracting the covariant decomposition (B.3) 
with 6364 and by using the relation 63 • 64 = — 1. The outcome is 

= -!liii±iiI^ + 0(e). (B.5) 

Integral lijlfi"^ {{q + pi) ■ 62)] 

The covariant decomposition of a rank-3 bubble reads as follows 



,.-li+Mli±|Iii±M: . {g^^p%+r^p^.+9^^P^:)Boo. 






with "pji = pj — Pi. After the contraction with g^'^e!^ we have 



PuP'nP'n Biiu (B.6) 



lij [^J■ {{q + Pi) ■ 62)] = j d q jj-jy = 2e (pj^ ■ 62) Booi 

- 9 

ITT 

= ^ iPJ^ ■ 62) (pI - 2mj - 4m2) + 0{e) . (B.7) 

Integral lijkifJ-*] 

We start from the decomposition 

j^. , I ,.^ i^ + P^nq + P^2q+J.nq + P.r ^ (^,.^p. ^ ^,p^.. ^ ^,.^.p) ^^^^^ 

+ (rank-4 tensors containing pji ,Pjk) (B.8) 
After the contraction with g^^" gP'^ the tensors containing pji, pjk vanishes and we get 

I^,k [^i'] = I d'q-^-^^ = 4e(e - 1) Coooo 

Pjk 'Pji ' Pki 2 2 2 I /-o/ ^ fo n\ 

^^ j m,, - m^ - rrif, + 0{e) . (B.9) 
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